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ABSTRACT 

Recurrent nova CI Aquilae is still bright 300 days after the optical maximum, showing the slowest 
evolution among recurrent novae. We predict the turn-off time of CI Aql 2000 outburst coming in August 
2001 after a supersoft X-ray source (SSS) phase lasts 250 days. We also predict helium enrichment 
of ejecta, He/H^ 0.25 by number. Observational confirmations are urgently required. Based on the 
optically thick wind mass-loss theory of the thermonuclear runaway model, we have also estimated the 
WD mass to be M W d = 1-2 ± 0.05 M© by fitting our theoretical light curves with the 1917 and 2000 
outbursts. The mass of the hydrogen-rich envelope on the WD is also estimated to be AM ~ 6 x 1O _6 M0 
at the optical maximum, indicating an average mass accretion rate of M acc ~ 0.7 x 1{)~ 7 Mq yr _1 
during the quiescent phase between the 1917 and 2000 outbursts. Using these obtained values, we have 
consistently reproduced the light curve in quiescence as well as the two outburst phases. Thus, we 
predict the turn-off time to be in August 2001 for the 2000 outburst. We strongly recommend soft X-ray 
observations to detect SSS until August 2001 because the massive wind phase have already ended in 
December 2000 followed by an SSS phase that very likely lasts until August 2001. 



Subject headings: binaries: close 
winds, outflows 

1. INTRODUCTION 



novae, cataclysmic variables — stars: individual (CI Aql) 
X-rays: stars 
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The second recorded outburst of CI Aquilae was dis- 
covered in April 2000 by Takamizawa (2000) since the 
first recorded outburst in 1917 (Reinmuth 1925; Duerbeck 
1987). Thus, CI Aql has been recognized as a member of 
the recurrent novae (Yamaoka, Ayani, & Shirakami 2000) . 
After optical brightness reached its maximum V ~ 9 in 
early May 2000, it rapidly decreased to V ~ 13 in about 50 
days. A plateau phase follows; the brightness levelled off 
at V ~ 13.5. CI Aql is still as bright as V ~ 14 until now 
about 300 days after maximum, being the longest record of 
on-time among the recurrent novae (e.g., VSNET archives, 
http : //www.kusastro .kyoto-u. ac . jp/vsnet/). 

Recurrent novae constitute a small class of cataclysmic 
variables, the outbursts of which repeat in a decade to 
a century. Ejecta of the recurrent novae are not metal- 
enriched but similar to the solar value [Z ~ 0.02), indicat- 
ing that the WDs in the recurrent novae are not eroded 
but have grown in mass (e.g., Starrfield, Sparks, & Tru- 
ran 1985). Masses of white dwarfs (WDs) in the recur- 
rent novae are estimated to be so massive as the Chan- 
drasekhar mass limit from the light curve fitting (e.g., 
Hachisu & Kato 1999; Hachisu et al. 2000a; Hachisu 
& Kato 2000a, b, 2001). These features strongly suggest 
that the recurrent novae are an immediate progenitor of 
Type la supernovae (SNe la) if the WD consists of carbon 
and oxygen (Nomoto, Thielemann, & Yokoi 1984). Evo- 
lutionary paths to SNe la via recurrent novae have been 



pointed out by Hachisu, Kato, & Nomoto (1999a). 

Recurrent novae are divided into three subclasses de- 
pending on their orbital periods: T Pyx subclass with a 
dwarf companion, U Sco subclass with a slightly evolved 
main-sequence (MS) or subgiant companion, and RS Oph 
subclass with a red giant companion (e.g., Warner 1995). 
CI Aql has an orbital period of 0.618 days, belonging to 
the U Sco subclass. It is therefore likely that the light 
curve of CI Aql develops like U Sco. 

The 1999 outburst of U Sco was densely observed from 
the rising phase to the final decay phase through the mid- 
plateau phase. The detection of supersoft X-rays during 
the plateau phase (Kahabka et al. 1999) strongly supports 
the theoretical description of nova outbursts — massive 
wind phase, steady hydrogen burning phase (no winds), 
and final cooling phase, — which corresponds to the rapid 
decline, plateau, and final decline phases in the light curve 
(e.g., Hachisu et al. 2000a), respectively. Because the 
development of CI Aql outbursts is much slower than that 
of U Sco, there is still a chance to observe supersoft X- 
rays. Supersoft X-rays are direct evidence of hydrogen 
shell burning on a WD so that they provide us valuable 
information on recurrent novae and SN la progenitors. 

In this Letter, we determine the mass and the growth 
rate of the WD as accurately as possible by theoretically 
analyzing light curves. Then, (1) we predict the turn- 
off time of nuclear burning and the duration of supersoft 
X-ray source (SSS) phase of the present outburst and (2) 
examine the possibility whether or not CI Aql will explode 
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as an SN la. In §2, we describe our model of CI Aql in 
quiescence and estimate the inclination angle, the size and 
thickness of the accretion disk, and so on. In §3, the 1917 
outburst is numerically analyzed and we determine the 
mass of the WD. In §4, the 2000 outburst is also numer- 
ically analyzed and the duration of wind is estimated. In 
§5, we predict the duration of the SSS phase. Discussion 
follows in §6. 

2. BASIC MODEL IN QUIESCENCE 

CI Aql is an eclipsing binary. We have considered a bi- 
nary model as illustrated in Figure 1 which consists of (1) 
the WD photosphere, (2) the MS photosphere which fills 
its Roche lobe, and (3) the accretion disk (ACDK). We 
assume that the size and thickness of the ACDK are given 
by two parameters, a and (3, as 

-Rdisk = aRl, (1) 

and 

h = (3R disk (-^-) , (2) 

V -Kdisk / 

where i?disk is the outer edge of the ACDK, Rl the effective 
radius of the inner critical Roche lobe for the WD compo- 
nent, and h the height of the surface from the equatorial 
plane, w the distance on the equatorial plane from the 
center of the WD. Here, we adopt tu-square law [y = 2) 
to mimic the effect of flaring-up at the rim of the ACDK 
(e.g., Schandl, Meyer-Hofmeister, & Meyer 1997). 

The accretion luminosity of the WD (e.g., Starrficld, 
Sparks, & Shaviv 1988), the viscous luminosity of the 
ACDK, and the irradiations by the WD (e.g., Schandl et 
al. 1997) are also included. Here, we assume that the 
surfaces of the MS and of the ACDK emit photons as a 
black-body at a local temperature of each surface patch 
area (see Fig. 1) which is heated by radiation from the 
WD. The efficiencies of irradiation of the MS and of the 
ACDK arc assumed to be 7/ms = 0.5 and t?dk = 0.5, re- 
spectively. The numerical method is essentially the same 
as that in Hachisu et al. (2000b) and fully described 
in Hachisu & Kato (2001). A circular orbit is assumed. 
Light curves of CI Aql in quiescence has been obtained by 
Mennickent & Honeycutt (1995). They determined the 
orbital period and ephemeris to be JD 2,448,412.167(25) 
+ 0.618355 xE at eclipse minima. 

Figure 2 shows the averaged (smoothed) observational 
points (open circles) together with our calculated V light 
curves for the suggested WD mass of Mwd = 1.2M© 
(i?WD = 0.0072i? Q ) and mass accretion rate of M acc ~ 
1 x 10~ 7 M Q yr _1 (see the following sections). To fit our 
theoretical light curves with the observational points, we 
have calculated V light curves by changing the parameters 
of a = 0.5—1.0 by 0.1 step, /? = 0.05—0.50 by 0.05 step, 
^ P h,MS = 5000—8000 K by 100 K step, T phjdisk = 3000— 
6000 K by 500 K step, and i = 70—90° by 1°' step and seek 
for the best-fit model by minimizing RMS of the residu- 
als. The light curves are also calculated for four compan- 
ion masses of Mms — 0.8, 1.1, 1.5 and 2.0M Q . Since we 
obtain similar light curves for all of these four cases, we 
adopt Mms = 1-5M0 because the mass transfer is as high 
as ~ 1 x 10~ 7 M Q yr _1 and this indicates a thermally un- 
stable mass transfer (van den Heuvel et al. 1992), which 
requires M M s/M WD > 1.1 (Webbink 1985). A 1.5 M MS 



expands to fill the Roche lobe of ~ 1.7 Rq after a large 
part of hydrogen at the center is consumed (e.g., Hachisu 
et al. 1999b). 

In the best-fit model, a rather blue color index of 
(B — V) c w 0.14 is obtained as shown in Figure 3. This 
suggests a large color excess of E(B — V) = (B — V) — 
(B - V) c w 0.86 with the observed color of (B - V) a ~ 1.0 
(Mennickent & Honeycutt 1995). This color excess is very 
consistent with E{B — V) = 0.85 estimated by Kiss et al. 
(2001). 

3. THE 1917 OUTBURST 

First we reproduce the light curve of the 1917 outburst. 
Its Tycho B magnitude light curve has recently been re- 
ported by Williams (2000), which is shown in Figure 4 
together with our calculated B light curves. Our numer- 
ical method for calculating light curves is fully described 
in Hachisu & Kato (2001). 

It has been established that the WD photosphere ex- 
pands to a giant size at the optical maximum in nova out- 
bursts. Since the separation of CI Aql is a ~ 4 Rq (Greiner 
& Alcala 1996), the WD envelope engulfs the MS compan- 
ion so that the early phase light curve is calculated only 
from the WD photosphere (e.g., Kato & Hachisu 1994; 
Kato 1999). The accretion disk and the companion do not 
contribute to the B light until it decreases to B ~ 13 mag. 
We determine the WD mass to be M WD = 1.2 ± 0.05 M Q 
from the early phase light curve as seen in Figure 4. 

Assuming Mwd = 1-2 M Q , we estimate the envelope 
mass at the optical maximum to be AM = 5.8 x 10 _6 M Q 
for the hydrogen content X = 0.7. If the duration of 
quiescence before the 1917 outburst is almost the same 
as the quiescent phase between the 1917 and 2000 out- 
burst, this envelope mass indicates a mass accretion rate 
of M acc ~ 0.7 x 1O~ 7 M yr" 1 if no WD matter has been 
dredged up. In our wind model, about 81% (4.7x 1O _6 M ) 
of the envelope mass has been blown in the wind while the 
residual 19% (1.1 x 1O~ 6 M ) has been left and added to 
the helium layer of the WD. The net mass-increasing rate 
of the WD is M Hc = 0.13 x 10~ 7 M Q yr" 1 . Thus, CI Aql 
may become an SN la if the donor is massive enough. 

4. THE 2000 OUTBURST 

CI Aql bursted in April 2000. Its optical maximum 
(V ~ 9 mag) was reached on May 5, 2000 (HJD 2451669.5, 
Kiss et al. 2001). As shown in Figure 5, the visual bright- 
ness quickly decreased to 13.5 mag in about 60 days. Then, 
it stays at V ~ 14 mag, i.e., a plateau phase. This mid- 
plateau phase is very similar to that of U Sco and can 
also be explained by the irradiation of the accretion disk 
(Hachisu ct al. 2000a). Therefore, we have reproduced the 
light curve in the plateau phase by assuming similar model 
parameters to those of U Sco. Theoretical light curves in 
Figure 5 are calculated for a pair of 1.2 M WD + 1.5 M 
MS. In the plateau phase of the light curve, i.e., when the 
WD photosphere becomes much smaller than the binary 
size, the light curve is determined mainly by the irradia- 
tions of the ACDK and of the MS. The unheated surface 
temperatures are fixed to be T p h,MS = 6300 K for the MS 
and to be T p h,di S k = 5300 K at the disk rim, which are the 
same as the basic model in quiescence. 
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The luminosity of the accretion disk depends on both 
the thickness (3 and the size a. We have examined a total 
of 160 cases for the set of (a, (i), which is the product of 
16 cases of a = 0.5 — 2.0 by 0.1 step and 10 cases of (3 = 
0.05 — 0.50 by 0.05 step. Here, we have adopted a set of 
a = 1.4 and /3 = 0.2 during the wind phase and a set of 
a = 0.7 and (3 = 0.2 after the wind stopped. 

The decline of the early phase (t ~ — 70 days) hardly 
depends on the hydrogen content X of the WD envelope 
but the ends of the wind phase and of the hydrogen shell- 
burning phase depend sensitively on the hydrogen con- 
tent X as shown in Figure 5. Matsumoto & Kato (2000) 
reported a sharp ~ 1 mag drop of R magnitude in late 
November 2000 and Kiyota (VSNET archives) also ob- 
served a ~ 1.5 mag drop of / magnitude as shown in Fig- 
ure 5. If we attribute this drop to the end of wind phase, 
the hydrogen content is somewhere between X = 0.35 and 
X = 0.5, more closer to X = 0.5. We here adopt X = 0.5. 
In the case of X = 0.5, the continuous wind stopped in 
early December 2000 about 220 days after maximum and 
the steady hydrogen shell-burning will end in early August 
2001 about 450 days after maximum as shown in Figure 
5. 

5. SUPERSOFT X-RAY SOURCE PHASE 

We here strongly encourage X-ray observations until the 
mid August 2001, that is, the end of steady hydrogen shell- 
burning. No supersoft X-rays are observed in the early 
decline phase because the photospheric radius of the WD 
is rather large and the photospheric temperature is rel- 
atively low. After the recurrent nova enters the plateau 
phase, the photospheric temperature of the WD increases 
to high enough to emit supersoft X-rays. During the mas- 
sive wind phase, however, we do not expect supersoft X- 
rays because they are self-absorbed by the wind itself. U 
Sco was observed as a luminous supersoft X-ray source in 
the mid-plateau phase of the 1999 outburst just after the 
massive wind stopped (Kahabka et al. 1999; Hachisu et 
al. 2000a). Our analysis suggests that the massive wind 
stopped in early December 2000 and that supersoft X-rays 
are now very likely detectable until August 2001. 

6. DISCUSSION 



The calculated distance modulus of the 1917 outburst is 
about (to — M)b = 14.5. Therefore, the distance to CI Aql 
is estimated to be 1.6 kpc for the reddening of E(B — V) = 

0. 85 and the absorption of A B = 3.5 and Ay = 2.7. The 
apparent distance modulus of (to — M)y — 13.7 from the 
2000 outburst as shown in Figure 5 is also consistent with 
the distance of 1.6 kpc and Ay = 2.7. 

On the other hand, we obtain an apparent distance mod- 
ulus of (m — M)v — 12.86 in the quiescent phase as shown 
in Figure 2, which corresponds to the distance of 1.1 kpc. 
This is not consistent with our estimation from the out- 
burst light curves. To examine this difference, we check 
other observations and found that three others indicate 

V « 16.2 in quiescence, which is 0.7 mag darker than 

V « 15.5 in Figure 1 of Mennickent & Honey cutt (see Ta- 
ble 2 of Mennickent & Honeycutt 1995; Szkody & Howell 
1992). The V magnitudes in their Figure 1 were obtained 
with a fully automated telescope and the errors in the zero 
point (formally 0.2 mag as in their text) could be larger 
than expected (private communication with Mennickent). 
If we adopt V = 16.2 in quiescence, then the distance 
modulus is (to — M)y = 13.56 and the distance to CI Aql 
becomes d — 1.5 kpc, being consistent with the distance 
estimate from the outbursts. Thus, we determine the dis- 
tance to CI Aql to be d ~ 1.6 kpc with Ay = 2.7. 

The lacks of flickering and of Balmer emission lines in 
quiescence lead Mennickent & Honeycutt (1995) to con- 
clude that CI Aql is not interacting at least at the present 
time. However, the inner part of the accretion disk cannot 
be seen from the Earth as easily understood from Figure 

1 . This is the reason for the lacks of emission lines and of 
flickering even when the system includes a mass-accreting 
WD. This is essentially the same situation as in U Sco 
(Schaefer 1990; Johnston & Kulkarni 1992; Schaefer & 
Ringwald 1995; Hachisu et al. 2000b). 

We are very grateful to the VSNET members who ob- 
served CI Aql and sent their valuable data to the VS- 
NET. We also thank R. Mennickent for his comment on 
their photometry and the anonymous referee for many 
kind comments that help us to improve the manuscript. 
This research has been supported in part by the Grant- 
in-Aid for Scientific Research (11640226) of the Japanese 
Ministry of Education, Science, Culture, and Sports. 
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Fig. 1. — Configuration of our CI Aql model in quiescence. The cool component (right) is a slightly evolved MS (1.5Mq) filling up its inner 
critical Roche lobe. Only north and south polar areas of the cool component are irradiated by the hot component (1.2 Mq WD, left). The 
hot component itself and the central part of the ACDK are not seen from the Earth because they are blocked by the flaring-up rim of the 
ACDK. Here the separation is a = 4.25-Rq; the effective radii of the inner critical Roche lobes are i?J = 1.53i?0, and ijj = R2 = l-69i?0, for 
the primary WD and the secondary MS, respectively. 
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1.2 M Q WD + 1.5 M MS T(MS) = 6300 K T(DK)=5500 K 
dM/dt=1.0E-7 a=0.7 0=0.3 

(m-M) v =12.86, 13.23, 13.43 for i=80°, 74°, 71° 
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Fig. 2. — Calculated V light curves are plotted against the binary phase (two phases from —1.0 to 1.0) together with the smoothed 
observational points (open circles) taken from Mennickent & Honcycutt (1995). The model is a binary system of 1.2Mq WD + 1.5Mq 
MS. The other parameters are printed in the figure. Light curves are plotted for different three inclination angles, i.e., i = 71° (dotted), 74° 
(dashed), and 80° (solid). 
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Fig. 3. — Calculated B — V light curve against the binary phase for the inclination angle of i = 80°. 
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Fig. 4. — Calculated B light curves are plotted against time (HJD 2,421,000+) together with the observational points of the 1917 outburst. 
Open circles indicate observational points (taken from Williams 2000). The model consists of a bloated WD photosphere with no accretion 
disk and a non-irradiated MS companion. The hydrogen content of the WD envelope is assumed to be X = 0.70 for all models. The apparent 
distance modulus of (m — M)b = 14.5 is assumed for all the WD masses. Solid lines indicate the light curves connecting the B light at 
the binary phase 0.5 (roughly the mean brightness) while dashed lines correspond to those connecting the B light at the binary phase 0.0 
(mid-eclipse). 
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Fig. 5. — Calculated V and B light curves are plotted against time (HJD 2,451,000+) together with the observational points of the 2000 
outburst. Small dots indicate observational V and visual magnitudes while open squares represent observational B magnitudes and open 
circles indicate observational I magnitudes (all taken from VSNET archives). Calculated V light curves are plotted for three different hydrogen 
contents of the WD envelope, X = 0.70 (thick dotted), X = 0.50 (thick solid), and X = 0.35 (thick dot-dashed). Thin solid line indicates 
£?-magnitude of X = 0.50 case. Each light curve connects the brightness at the binary phase 0.5 (roughly the mean brightness). The apparent 
distance moduli of (rre — M)y = 13.7 and (m — M)g = 14.5 are obtained by fitting. 



